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ABSTRACT 

Chandra j4C/S'-S3 observations of the nearby SO galaxy NGC 1332 resolve much of the X-ray emission 
into 73 point-sources, of which 37 lie within the D25 isophote. The remaining galaxy emission comprises 
hot, diffuse gas and unresolved sources and is discussed in two companion papers. The point-source 
luminosity function (XLF) shows the characteristic break seen in other early-type galaxies at ~ 2 x 10'^^ 
erg s""'^. After applying corrections for detection incompleteness at low luminosities due to source 
confusion and contamination from diffuse galactic emission, the break vanishes and the data are well- 
described as a single power law. This result casts further doubt on there being a "universal" XLF 
break in early-type galaxies marking the division between neutron-star and black-hole systems. The 
logarithmic slope of the differential XLF (dN/dL), (3 = 2.7±0.5, is marginally (~ 2.5(t) steeper 
than has been found for analogous completeness-corrected fits of other early-type galaxies but closely 
matches the behaviour seen at high luminosities in these systems. Two of the sources within D25 are 
Ultra-luminous X-ray sources (ULX), although neither have Lx> 2 x lO'^^ erg s""'^. The absence of 
very luminous ULX in early-type galaxies suggests a break in the XLF slope at ~l-2xl0'^^ erg s"'^, 
although the data were not of sufficient quality to constrain such a feature in NGC 1332. The sources 
have a spatial distribution consistent with the optical light and display a range of characteristics 
that are consistent with an LMXB population. The general spectral characteristics of the individual 
sources, as well as the composite source spectra, are in good agreement with observations of other 
early-type galaxies, although a small number of highly- absorbed sources are seen. Two sources have 
very soft spectra, two show strong variability, indicating compact binary nature and one source shows 
evidence of an extended radial profile. We do not detect a central source in NGC 1332, but we find 
a faint (Lx= 2 ± 1 x 10^^ erg s~^) point-source coincident with the centre of the companion dwarf 
galaxy NGC 1331. 

Subject headings: galaxies: elliptical and lenticular, cD — galaxies: individual (NGC 1332) — X-rays: 
binaries — X-rays: galaxies 



1. INTRODUCTION 

Prior to the launch of Chandra only a small num- 
ber of the very brightest point-sources in early type 
galaxies could be resolv e d from the diffuse g alactic 
emission (fF abbiano 1989; Colb ert fc Mushotzkvl 119991: 
iRoberts fc Warwick 2000). The study of their point- 
source populations was therefore restricted to the aver- 
age, composite properties inferred by decomposing the 
emission into a number of spectral or spatial model com- 
pon ents fe.g. iMatsushi ta et al.llT99^ iBrown fc BregmanI 
l2001ip . The advent of Chandra has revolutionized this 
field, however, allowing a large fraction of the sources 
to be re solved and studied directly (for a recent re- 
view, see lFabbiano fc White 2003). Extragalactic point- 
sources most probably represent an heterogeneous mix- 
ture of different source types, although the old stellar 
populations characteristic of early-type galaxies suggest 
a predominantly low-mass X-ray binary (LMXB) na- 
ture. As endpoints of stellar evolution, the properties 
of X-ray binaries are a crucial diagnostic for the evo- 
lution of the stars as a whole within the galaxy. Al- 
though the data are generally of insufficient quality to 
investigate fully for each source t he rich, diagnostic phe- 
nomenology of Galactic LMXB (|White et al.lll99,'iD . in 
some cases it has been possible to identify variability (e.g. 
I^razin et al. 2001; Kraft et al. 2001^ and the spectral 
signatures of black-hole binaries (e.g. iMakishima et all 



l2(i00tlHumnhrev et al.ll200l . 

In order to study the properties of the X-ray 
point-source population as a whole, it is common 
practice to consider the X-ray luminosity function 
fXLF^ (e.g. ISarazin et^aU HM IBlanton et a,lJ 
iZezas fc Fabbianol |2002() . the shape of which is a 
strong fun ction of the ag e of the stellar population 
l|Kilgard et al. 2002: Bel czvnski et al.ll2003D . Between 
early-type galaxie s, there is remarkable sim ilarity in the 
shape of the XLF ijKim fc Fabbianol2003b() . which is typ- 
ically found to be a steep power law with a break oc- 
curring around 2-4 xlO'^* erg s~^ and a high-luminosity 
slope, 2± 2-3 (.Sarazin et al. 20 01: Blanton et al. 200j 
iKraft et al.ll200lUColbert et aljl200i^ . These slopes are 
also br oadly cons istent with old stellar populations in 
M31 (Ko ng et al.l l2003'). It has been suggested that the 
presence of breaks at luminosities at around the Ed- 
dington limit of a 1.4Mq neutron star (Ledd =2-4x10"^^ 
erg s~^, depending on the composition of th e accreting 
matt er and the neutron-star equation of state: iPaczvhskil 
IT983I) may arise from a di vision between neutr on-star and 
black-hole systems (e.g. iSarazin et aD 12001 . However 
in the elliptical galaxy NGC 720, a substantially higher 
luminosi ty break, at ^ 1 x 1 ^^ erg s~^ was found in- 
stead bv lJeltema et al.ll|2003D . iKim fc Fabbianol l)2003afl 
pointed out that in the elliptical galaxy NGC 1316, a 
break is arises only if no correction is made for point- 



2 



Humphrey and Buote 



source detection incompleteness. App lying this correc- 
tion to a sample of early- type galaxies. iKim fc Fahhiand 
()20031:)j) found a similar result, although they still found 
marginal evidence of a break at ~ 5 x 10'^* erg s'-^. This 
luminosity is sufRciently high, however, to cast further 
doubt upon the idea that it marks the division between 
neutron-star and black-hole binaries. 

At the extreme end of the XLF are the so-called 
"Ultra Luminous X-ray sources" (ULX), which are 
non-nuclear objects with luminosities exceeding 10'^^ 
erg s ~^, somet imes reaching as high as 10"^^ erg s~^ 
(e.g.lFabbianolfl OSQ: Makishima et alJl2000t IZezas et all 
120021 IDavis fcMushotzkv..2003,l . They have been found 
in ga laxies of all morphological types ijColbert &: PtakI 
12002^. although there is a strong association between 
ULX and star-formation so that their sp ecific frequency 
is far lower in early-type galaxi es ( Kilg ard et alJ 120021: 
iHumphrev et al.l 120031: llrwin et 'al...2003bj). However, in 
the ot herwise normal elliptical NGC 720, IJeltema et al.l 
(IPOS') found a remarkable population of 9 ULXes. The 
nature of ULX is somewhat enigmatic. Although the less 
extreme objects (Lx~l-2xl0^® erg s~^) are consistent 
with isotropic emission from Eddington-limited black- 
hole binaries with compact-object masses ~1O-2OM0, 
there has been much debate over whether the most 
luminous objects represent "Intermediate mass black 
holes" with mases >1OOM0, super-Eddington emis- 
sion from sources in the thermal-timescale mass-transfer 
phase, or beamed emiss i on from lower-luminosity s ources 
ilColbert & Mushotzkvl IToOfll: IKing et all 120011 iKmel 
l200ir 

Another feature of the LMXB population of particular 
interest is its possible association with globular clusters 
(GCs), since it provides valuable insight into X-ray bi- 
nary formation. The fraction of the extragalactic LMXB 
identified with GCs varies consi derably between galaxies , 
from ^20-70% llAngelini et alJlMlt fKundu et alJl2002t 
iSarazin et al.ll200H) . Nonetheless, there is evidence that 
the frequency of at least part of the LMXB population 
correlates with the numbers of GCs so that ~4% of all 
GCs contain a n LMXB, consistent with the Milky Way 
ijKundu et allOoOSHSarazin et al.ll2003j) . These authors 
also found evidence that LMXB favour brighter, redder 
GCs. 

In this paper, we analyse with Chandra the hith- 
erto unresolved point-source population of the lenticu- 
lar gal axy NGC 1332. Usin g Rosat PSPC data of this 
galaxy, IBuote fc Canizare^ (^1996, hereafter BC) found 
elongation of the X-ray isophotes and demonstrated that 
the galaxy mass profile was marginally inconsistent with 
the optical light, indicating the presence of a substan- 
tial dark matter halo. Unfortunately, most of the point- 
source population could not be resolved, so that it was 
necessary to estimate limits for the unresolved point- 
source contribution to the X-ray emission, and tight con- 
straints on the shape of the da rk matter halo could not be 
obtained. Using ASCA data IBuote fc CanizaresI 1119971) 
found a hard spectral component in the integrated stellar 
emission, indicating a significant contribution from unre- 
solved LMXB, although this did not allow these sources 
to be studied in detail. Chandra allows us for the first time 
to separate a large fraction of the LMXB from the dif- 
fuse emission of NGC 1332, and to study them directly. 
The properties of the diffuse gas and the mass distribu- 



tion are discussed separately in two com panion papers 
l|Humphrev et all 120041: IBuote et al.ll2004 hereafter Pa- 
pers I and II, respectively). 

In order to determine the properties of the stellar 
population accurately, it is important to adopt a re- 
liable distance estimate (c.f. NGC 4038/9, for which a 
recently revised distance estimate would reduce the re - 
ported ULX population threefold; iSaviane et al.ll20d3j) . 
For NGC 1332, distances determined both by surface 
brightness fluctuations (SBF) and the globular cluster lu- 
minosity function (GCLF) are in excellent ag reement, be- 
ing 23 ± 3 Mpc and 22 zfc 6 Mpc, r espectively ijTonrv et alJ 
l2?)0lllKundu & Whitmoreil2001D . All errors quoted here, 
and subsequently, are 90% confidence limits unless oth- 
erwise stated. 

2. OBSERVATIONS AND DATA ANALYSIS 

The region of sky containing NGC 1332 was observed 
with the ACIS instrument aboard Chandra between 2002 
September 19 10:39 and September 20 02:59 UTC, the 
galaxy being centred on the S3 chip, for a nominal 
~60 ks exposure. An additional observation made be- 
tween September 18 02:56 and 08:27 UTC, for a nom- 
inal 20 ks exposure, was heavily contaminated by flar- 
ing of the background and so we present results only 
for the better-quality dataset. For data- reduction, we 
used the CIAO 3.0.1 and Heasoft 5.2 suites of soft- 
ware tools, and for spectral analysis we used XSPEC 
11.2.0 and ISIS 1.1.5. In order to use the most up-to- 
date calibration, we reprocessed the Level-1 data, fol- 
lowing the standard Chandra threads^. We corrected a 
time-dependent drift in the gain with the apply_gain^ 
task. To identify periods of unusually high background 
for which the degradation in signal-to-noise can affect 
analysis, a lightcurve was generated from source-free re- 
gions of the active chips and examined by eye for "flares" . 
For much of the observation there was evidence of low- 
level flaring (AR/R ~20%, where R is the count-rate). 
For point-source identification and characterisation, we 
found the S/N improvement from reducing the back- 
ground level (by ~10%) was outweighed by the reduc- 
tion in the exposure-time (by ^30%) if we excised data 
during the strongest flares. We therefore used the entire 
observation, giving a total exposure time of 57 ks. We 
have systematically confirmed that including the flaring 
did not bias our results. A more thorough discussion of 
the flaring is given in Paper II. Fig.^shows the image of 
the inner part of the S3 chip, which has been flat-fielded 
with an exposure-map generated with CIAO at 1.0 keV 
(approximately the peak energy of the diffuse emission) , 
and smoothed with the CIAO task csmooth (taking the 
default parameter values). A number of bright point- 
sources are clearly visible in the X-ray data. For com- 
parison purposes the Digitized Sky Survey (DSS) optical 
image is also shown. 

3. POINT SOURCE DETECTION 

Point-source detection was performed using the CIA 
tool wavdetect (Freeman et al. 2002). In order to im- 
prove the likelihood of identifying sources with peculiarly 
hard or soft spectra, full-resolution images were created 

^ http://cxc.harvard.edu/ciao/threads/index.htmll 
^ httpT/cxcTarvard^au/con 
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Fig. 1. — left panel: False colour, smoothed image of the central part of the S 3 chip. The image was initi ally flat-ficldcd with an 
exposure- map created at 1 keV (peak energy). Overlaid is the D25 ellipse taken from lde Vaucouleurs et al.l 119911) . but the position angle 
has be en altered to 116° to agree with BC. Right panel: the DSS optical image, overlaid with the D25 ellipse and, following Humohrcv ct sdj 
120031) . the positions of X-ray sources with Lx> 10^^ erg s~^ (X) and those with 5 X 10^* <Lx< 10^^ erg s~^ The identification 

numbers of a selection of interesting sources are given on the figure. The bright galaxy to the left of NGC 1332 is NGC 1331. 



of the region of the ACIS focal-plane containing the S3 
chip in three different energy-bands (0.1-10.0 keV, 0.1- 
3.0 keV and 3.0-10.0 keV). Sources were detected sep- 
arately in each image. In order to minimize spurious 
detections at node or chip boundaries we supplied the 
detection algorithm with exposure-maps generated at en- 
ergies 1.7 keV, 1.0 keV and 7 keV respectively (although 
the precise energies chosen made little difference to the 
results). The detection algorithm searched for structure 
over pixel-scales of 1, 2, 4, 8 and 16 pixels, and the de- 
tection threshold was set to ^ 10^^ spurious sources per 
pixel (corresponding to ^0.1 spurious detections per im- 
age). The source- lists obtained within each energy-band 
were combined and duplicated sources removed, and the 
final list was checked by visual inspection of the images. 

A total of 76 sources were detected on the S3 chip, 
of which the brightest and most extended was cen- 
tred within 1" of the galaxy centroid as given in NED 
(which is accurate to ^ 2"). Examination of the 
azimuthally-averaged radial brightness profile of this 
"source" (Sect. 14 .Sf) clearly indicated it to be very ex- 
tended and so we identified this with the central part 
of the diffuse galactic emission (which had been spu- 
riously identified by wavdetect as a point-source due 
to its being very centrally-peaked). We found that one 
point-source lay within the B-band 25*^^ magnitude (02.=^) 
isophote of NGC 1331 (as hsted in lde Vaucouleurs et al.l 
^91). In fact, this source was within 1" of the 
galaxy centroid listed in NED, but was rather faint 
(with a count-rate ~ 4 x lO^^count s~^, corresponding 
to Lx ~ 2 ± 1 X 10'^* erg s~^). In addition, one source 
had no photons in the 0.5-7.0 keV band, suggesting that 
it was spurious and we therefore omitted it from our 
list. This left 73 sources on the S3 chip, of which 37 
lay within the D2.'s i s ophot e of NGC 1332 as given in 
Ide Vaucouleurs et al.l ljl99H) . which has semi- major and 
semi-minor axes of 2.3' and 0.7', respectively. We ad- 
j usted the position angle given in |de Vaucouleur s et aP 
to 116° to agree with BC and match better the 



optical data. The source positions are given in Tabled 
It is interesting to compare our detections with those 
made with Rosat. Using the PSPC, 5 point sources were 
resolved within the central ring (BC). There are bright 
Chandra detections (which we labelled sources 1, 2 and 
6 in Table lU within 15" (- the PSPC point-spread) of 
three of these sources. One of the Rosat sources would 
not fit on the S3 chip and the remaining source was close 
to only one Chandra detection, source 71, which was too 
faint (Lx^ 10'^* erg s~^) to have been seen with Rosat. 
This indicates that the source seen with Rosat (which 
may, or may not, correspond with this Chandra detec- 
tion) exhibited dramatic variability over the ~ 6-year 
period between the observations. 

For each source, lightcurves, spectra and spectral re- 
sponses were generated using standard CI AO tasks. Lo- 
cal background regions were chosen for each source to 
ensure that contamination from the diffuse galaxy emis- 
sion (or from the mild flaring) did not bias our re- 
sults. These were accumulated from annuli centred on 
the source which were truncated so as to be contained 
entirely within the same ACIS node. They were cho- 
sen to cover an area at least 8 times the extraction re- 
gion, and containing at least 50 photons. We excluded 
from our background computation photons within a re- 
gion around each source obtained by scaling the axes of 
the appropriate 1-cr encircled-energy ellipse by a factor 6 
(we excluded such a large region since a few percent of 
a source's photons will leak from a region even as large 
as the 3-(T encircled-energy ellipse, which is problematic 
in the vicinity of bright sources). To compensate for the 
progressive quantum-efficiency degradation at low ener- 
gies due to contamination on the optical blocking filter, 
we used the apply_acisabs'^ script to correct the an- 
cillary response files. In order to investigate the possi- 
ble extent of these sources, dedicated software was used 
to generate radial brightness profiles centred upon each 

^ see|http://cxc.harvard.edu/ciao/threads/apply_acisabs| 
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Fig. 2. — Hardness ratios for all sources on the S3 chip with more 
than 40 counts. Overlaid is a grid of the loci of a simple absorbed 
power law model as T and Njj vary. F increases upwards, with grid 
lines shown for F=0. 0-3.0 in steps of 0.5. Several representative 
lines are labelled. Njj increases to the left and grid-lines are shown 
for Nh=0.0, 0.25, 0.5, 1.0, 2.0, 4.0 and 8.0 xlO^^cm'^, where Nh 
is column in excess of the appropriate Galactic value. Several lines 
of constant Njj are labelled on the figure in units of lO^^cm"^. 
The identification numbers of a selection of interesting sources are 
indicated on the figure. Sources inside the D25 region are shown 
as triangles, and those outside it are shown as circles. 



source and excluding photons from the vicinity of the 
other detected point-sources. The profiles were extracted 
initially in very fine, equally-spaced bins (^--^0.5" width), 
out to radii as large as 1' in order that a reliable back- 
ground could be estimated. For those sources close to the 
centre of the galaxy we adopted smaller radii (down to 
~0.5') to minimize contamination from the diffuse gas, or 
subsequently included a model to account for the diffuse 
component when fitting the profile fSect 14. 5]! . 

4. POINT SOURCE PROPERTIES 
4.1. Hardness ratios 

In order to compare the spectral phenomenology of in- 
dividual low-count sources it is convenient to construct 
hardness ratio s (e.g. White & MarshaU 1984: Kim et aj 
ezas et al. 2002). Background-subtracted counts 
were measured for each source in three energy-bands, a 
soft-band (S) from 0.3-1.5 keV, a medium-band (M) from 
1.5-3.0 keV and a hard-band (H) from 3.0-5.0 keV, and 
the hardness ratios HRl and HR2 were computed. These 
are defined, respectively, as HRl = (S — M)/(S + M) and 
HR2 = (M - H)/(M H). With these definitions, HRl 
is principally sensitive to the hydrogen column-density in 
the line-of-sight, and HR2 depends upon the steepness of 
the spectrum, allowing us to construct a useful "colour- 
colour" diagram (i.e. HR2 versus HRl) for NGC 1332 
(Fig|2J). In order to show how the hardness ratios de- 
pend upon the underlying spectrum, we have overlaid a 
grid showing the loci of simple absorbed power law spec- 
tra with a variety of different value of Nh and T. The 
position of a given source relative to the grid has implica- 
tions for its nature; most neutron-star or low-state black 
hole binaries, for example, would be expected to lie be- 
tween the r = 1 and the F = 2 lines, whereas high state 
black holes may be exp ected to lie above the F = 2.5 line 
l)Humphrev et al.l l200m. 

In NGC 1332, most of the sources are found 
to occupy regions consistent with Nh— Galactic 
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Fig. 3. — Composite spectrum for all sources in the D25 region, 
folded through the instrumental response, and shown with the best- 
fit bremsstrahlung model. 



ijPickev fc LockmanllT990|) and F ~l-2, as seen in other 
early- type galaxies. There is considerable scatter, but 
this merely refiects the few photons detected in each 
source. Several sources show some evidence of intrin- 
sic, or local, absorption in excess of the Galactic column 
but there are only two clear outliers (Sources 4 and 21). 
These are confirmed by spectral-fitting; e.g. Source 4 has 
Nh— 10^^ cm~^, in agreement with its hardness-plane 
position. Sources 6 and 21 have rather soft spectra, al- 
though both sources are somewhat harder than the "Su- 
persoft" and "Quasisoft" sources seen in other galaxies 
ijDiStefano fc Koii3l2003l) . Based on their colour-colour 
plane loci, they may be consistent with high-state black 
hole systems, although this is only one possibility, which 
will need independent confirmation. Individual spectral- 
fitting of these sources does not allow the spectra to be 
constrained very well, but if free, the F of a power law 
fit to each source spectrum does, indeed, rise to >3.0, 
consistent with this picture. 

4.2. Composite spectra 

In order to examine the overall properties of the source 
population, and specifically to allow us to separate the 
unresolved source component from the hot gas during 
spectral analysis of the diffuse emission (Paper II) , we ac- 
cumulated a composite background-subtracted spectrum 
of all the point-sources. Response matrices were gener- 
ated for the data using the mkwarf and mkrmf CI AO 
tasks, and corrected for the effects of quantum efficiency 
degradation at low energies analogously to the individual 
source spectra. Since all of the brightest sources detected 
in NGC 1332 were outside the D25 region, we separately 
constructed a composite spectrum for only those sources 
within this region (thereby preventing our results from 
being unduly biased by the brightest objects). We re- 
binned the data to ensure that in each bin there were 
at least 20 photons (to allow the use of the statistic) 
and a signal-to-noise ratio exceeding 3. Both spectra 
could be fitted excellently with simple one-component 
bremsstrahlung or power law models, with Nh fixed at 
the nominal Galactic value, as shown in Table |21 Only 
for the power law fit to all the S3 sources did freeing the 
Nh led to a modest improvement in the fit statistic. The 
spectra of the sources within the D25 region and those on 
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the entire chip were only slightly different. For our fit of 
sources within the D25 regio n, our results were in excel- 
lent agreement with those of llrwin et a l] ll2003ar . consis- 
tent with the apparently "universal" nature of the com- 
posite spectra of the lower-luminosity sources in early- 
type galaxies. 

4.3. Fluxing 

To measure the fluxes of the sources, we adopted two 
complementary techniques (which are effectively equiva- 
lent for the lower-counts sources) . We rebinned the spec- 
trum of each source to ensure at least 20 photons in each 
spectral bin (thereby allowing the use of x^)- When this 
gave fewer than 4 spectral bins, we instead regrouped 
the spectrum into ~12 bins of approximately equal width 
in energy and adopted the Cash (maximum likelihood) 
fit statistic (determining the goodness-of-fit from Monte- 
Carlo simulations). Initially we fitted a simple power law 
model, the slope of which matched the co-added source 
spectra within the D25 region. When the "null hypothe- 
sis probability" (i.e. the probability of the model fitting 
the data) was only acceptable at the 5% level or less, 
we rejected the model and allowed, separately, the pho- 
ton index and the Nh to be free. In all cases we were 
then able to obtain statistically acceptable fits. Although 
more complex spectral mo dels ar e generally more appro- 
priat e for Galactic LMXB l|Churc h & Bahicihska- ChurchI 
1200 ID . a simple power law can be used to provide an accu- 
rate flux estimate in a relatively narrow energy-band and, 
for so few photons, more complicated models could not 
be constrained. To determine the flux errors, we mod- 
ified the existing XSPEC-models to redefine normaliza- 
tion as the flux within the 0.3-7.0 keV band. With this 
deflnition, flux errors could be obtained trivially without 
correlations between different flt parameters biasing our 
results. For most sources the data were consistent with 
the co-added spectrum, although several of the brighter 
objects required additional Njj and in two cases we found 
an appreciably softer spectrum. Whilst this method pref- 
erentially corrects the brightest sources, only in a few 
cases was there a dramatic change in the Lx- Consid- 
ering the large luminosity errors of the faintest sources, 
this is unlikely to bias our results substantially. 

In order to compare more directly with other au- 
thors, we additionally computed counts-to-flux conver- 
sion factors for each of our sources. To allow for 
effective-area variations between detections, we folded 
the best-flt composite source hremsstrahlung model 
(NH=0.02xl022cni-2^ kT=8 keV) through the instru- 
mental responses of each source to determine a local 
counts-to-flux conversion factor. The Poissonian errors 
upon the measured counts were estimated using the 
Gehrels' approximation and translated into flux errors. 
We found excellent agreement between both methods of 
flux estimation, except in a few sources which showed 
unusually absorbed spectra. The results of our flux es- 
timation are shown in Table ^ which also deflnes our 
source naming convention (sources being numbered in 
descending order of luminosity) . 

For the distance to NGC 1332 (23 Mpc), we found that 
8 sources have Lx> 10^^ erg s~^, which qualifles them as 
ULX. The expected number of bright background objects 
in this sample is, however, appreciable (4, based on the 
Chandra deep fleld observations; Sect. 0), so some are 
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Fig. 4. — Full energy-band "lightcurve" of the variable source 51. 
Since so few photons were detected, we adopt a non-standard way 
of showing the lightcurve. In the lower panel, each arriving pho- 
ton is flagged as a vertical line, whereas the upper panel shows a 
comparison between the measured cumulative arrival probability 
of the photons in the source (solid line) and that expected for a 
non-varying source (dashed line) 



undoubtedly background objects. In order to compare 
with other authors, we considered separately only those 
ULX within the D25 region, of which there arc 2, with 
an expected background number of 0.3. Neither of these 
sources has Lx> 2 x 10^^ erg s~^. 

4.4. Variability 

To search for variability we accumulated source and 
background lightcurves with the maximum temporal res- 
olution (3.24 s) for each source. Since the time bin- 
ning was very flne (Tbin ^ Texp), we were able to com- 
pare the observed distribution of photon arrival times 
with that expected for a non-varying source, using 
the Kolmogorov-Smirnov (K-S) test (which requires un- 
binned data). The K-S test is fairly sensitive to large- 
amplitude variability occurring over timescales compa- 
rable to Texp. To minimize false positive detections, 
since we were considering a total of 74 sources (including 
the point-source in NGC 1331), we required that the null 
hypothesis (no variability) be rejected at 99.93% signifl- 
cance (corresponding to a 5% probability of a single false 
positive). 

We found that, of our sources, only two were variable at 
the required significance (as indicated in Table Since 
the background is known not to be entirely constant (due 
to fiaring), we tested these lightcurves having excluded 
flaring and found a comparably strong detection of vari- 
ability. In addition, we further tested the entire locally- 
accumulated background lightcurves of these sources us- 
ing the same test. In neither case could the null hypoth- 
esis (no variability in the background) be rejected at bet- 
ter than 7% signiflcance. It is therefore highly improba- 
ble that the measured variability of the sources arises 
from background fluctuations. Considering only the 
sources within the D25 region (and consequently lowering 
our detection signiflcance threshold) did not lead to any 
further detections of variability. Both variable sources 
had luminosities ~ 1-3 xlO'^* erg s~^, comparable to the 
Galactic "Z-track" sources, which can show variability 
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over similar time-scales (jHasinger fc van der Klialll98l . 
There were too few photons from each source, however, 
to classify these objects definitively. This will require a 
more sensitive, long-term study to track the spectral and 
timing behaviour. 

4.5. Spatial extent and the search for a central source 

The detection algorithm is capable of identifying both 
point-like and extended sources. In order to identify any 
sources which may be extended, we examined the radial 
brightness profiles accumulated for each source. Fitting 
was performed using dedicated software which was able 
to make use of the radial profiles described in Sect.|31and 
to take account of the instrumental point-spread at the 
peak photon energy of each source. The source data were 
rebinned into ~ 90 logarithmically-spaced radial bins and 
then adjacent bins were regrouped to ensure none con- 
tained fewer than 20 photons. The data were initially fit- 
ted with a model comprising a constant background term 
and a point-source at the extraction centroid. To account 
for exposure variation over the S3 chip, an exposure-map 
was generated for the source image at 1.7 keV and used 
to correct the applied model. For those sources within 
~0.5' of the galaxy centroid, the substantial diffuse emis- 
sion contaminated the profile sufficiently that a constant 
background model was no longer tenable. In these cases, 
we included an additional component, corresponding to 
a /3-model (with parameters matching the diffuse galac- 
tic emission; Paper III) and centred at an appropriate 
distance from the point-source. 

In all except one case we found no evidence of spatial 
extent greater than the instrumental PSF, placing limits 
of ^^50-100 pc (depending on the position in the field-of- 
view) on the actual size of these objects. Fig.^lshows the 
results for Source 14, which is clearly more extended than 
the PSF. Replacing the point-source with a /3-model im- 
proved the x^/dof from 104/60 to 59/58, which is signifi- 
cant at > 99.999% (upon the basis of an f-test). We were 
also able to fit the radial profile as two point-sources, 
one centred at the extraction centroid and one l.O^Q gg" 
away from it, of approximately equal brightness. Since 
we might expect two sources this far apart to be dis- 
tinguishable by Chandra, we tested this conclusion by 
simulating an image of two-point sources (separated by 
1") at an equivalent position in the focal-plane (gener- 
ating 1.0 keV point-source images from the Caldb PSF 
hypercubes). If there were more than ~ 100 counts in 
the sources they could be resolved. When we reduced the 
number of photons to ^ 25 each (corresponding to the 
observed number of counts in Source 14), Poisson noise 
blurred the outlines so that they remained unresolved. 
This is a practical example of source-confusion. The 
probability of a chance alignment of a background object 
within 1" of one of our detected sources was, however, 
only ~1%. Such a low probability does tend to suggest 
that the source was genuinely extended (having a size ~ 
1"), although we cannot altogether rule out the possibil- 
ity of a chance superposition. Future observations may 
be able to resolve the issue by identifying variability in 
this source, which would place limits on its size. 

In addition, we searched for a point-source embed- 
ded in the diffuse emission at the galactic centre. The 
radial profile of the diffuse emission out to 4.2' was 
grouped into ^150 logarithmically-spaced bins. These 



were regrouped to ensure at least 20 photons in each 
bin. The data could be well-approximated with a sin- 
gle /3-model plus a constant background term (folding 
in the instrumental PSF at 1.0 keV and correcting the 
model for exposure variations and gaps arising from re- 
moved point-sources), as shown in Fig. |31 Although the 
profile is extremely peaked (Rcorc = 1-2 ± 0.1") its cen- 
tre is significantly broader than the PSF of ACIS at 
its position. Adding a central point-source model does 
not improve the fit quality, placing an upper limit of 
13.4x10"'' count s"^ (Lx < 8.4 x 10^^ erg s~^) on any 
central source. A more detailed discussion of the mor- 
phology and radial profile of the diffuse emission will be 
given in Paper III. 

4.6. Source identification 

In order to identify possible counterparts to the point- 
sources we checked for objects listed in NED, Simbad 
and the Tycho-2 catalogue of Galactic (foreground) stars 
ffiag et al.ll200g i within 2" of each source. We did not 
identify any obvious counterparts to X-ray sources which 
would enable us independently to check the absolute as- 
trometry. However, the centroid of the diffuse galac- 
tic emission of NGC 1332 was located within 1" of the 
2MASS galaxy centroid (listed in NED), which is accu- 
rate to ~2" and there are no known pointing problems 
reported for our Chandra data. Therefore, we assumed 
that the astrometry of the detected sources was accurate 
to 2". Only one counterpart was found to any source, 
that being NGC 1331, the centroid of which lies within 
1" of an X-ray source (which we have not listed in Ta- 
ble [TJ. 

In addition, we checked for coincidence between the 
point source positions and a list of GC candidates kindly 
supplied by A. Kundu (priv. comm.). The GC candi- 
dates were identifi ed in HST WFPC2 observ ations cen- 
tred on NGC 1332 UKundu & Whitmord20niT) . which en- 
closes only ^50% of the D25 region. In order to perform 
this comparison it was necessary to check the relative 
astrometry between Chandra and HST. A comparison of 
the X-ray centroid and the optical centre in the appropri- 
ate archival WFPC2 image clearly revealed a systematic 
offset of '^l.G" between Chandra and HST. Compensat- 
ing for this effect, we found that 9 sources lay within 
0.75" of a GC candidate and so we assumed they were 
associated with a globular cluster. We adopted this limit 
since it corresponded to an apparent break in the distri- 
bution of offsets between GC and the X-ray source po- 
sitions. A total of 30 X-ray sources coincided with the 
WFPC field-of-view (of which 23 were within D25), and 
there were a total of 204 GC candidates in this region 
of NGC 1332, so that ~30% of the X-ray sources within 
this region were identified with GCs and 4.4% of all GCs 
hosted LMXB. Both of these results are c onsistent with 
observatio ns of oth er early-type galaxies ijKundu et alJ 
2002; Sara zin et aljf 2003. e.g.). We investigated whether 
there was any evidence of a systematic difference in the 
luminosities of the GC and the non-GC sources by com- 
paring the fluxes of the two populations using a two- 
sample Kolmogorov-Smirnov test. We found no evidence 
of a statistically significant discrepancy between the two 
distributions, which agreed with ^60% probability. The 
combined spectra of the GC sources exhibited no evi- 
dence of a systematic difference from the other sources. 
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Fig. 5. — Left panel: Full energy-band radial brightness profile for source 14. Three fits to the radial brightness profile are also shown. 
These comprise a constant background term and, respectively, the expected satellite point-spread computed at the peak source energy 
of 1.3 keV (solid line), a beta-model (dash-dotted line) and a pair of point-sources separated by 1" (dashed line). Right panel: radial 
brightness profile of the diffuse central emission from 0.3-7.0 keV. The data are shown with a fit comprising a constant background term 
plus a simple /3-model. The /3-model component is also shown separately (dash-dot line), and, for comparison, the azimuthally-averaged 
PSF (dotted line). There is no evidence of a central point-source. 



4.7. Spatial distribution of sources 

We examined the spatial distribution of the point- 
sources in order to determine to what extent they may 
follow either the optical light or the X-ray surface- 
brightness of the gas. In order to do this, we initially 
accumulated a histogram of the number of sources within 
each of 128 evenly-spaced radial bins centred at the 
galaxy centroid, and excluding the inner ~ 0.2' (where 
source detection incompleteness is most severe). In or- 
der to improve the signal-to-noise level, the data were re- 
binned until each bin contained at least 5 sources. Mod- 
els were tested against the data using the Cash-statistic 
(the goodness-of-fit being determined by Monte-Carlo 
simulation). We also tried using the statistic (which 
was less rigorous due to the low count-statistics) and ob- 
tained comparable results. We tested two different mod- 
els against the data, one following the optical light of the 
galaxy and one following the X-ray contours. To approx- 
imate the optical profile it was convenient to fit a single 
de Vaucouleurs profile with its effective radius fixed to 
31". We chose this value by correcting the semi-major 
axis effective radius measured by BC to account for the 
fact that our data were azimuthally averaged and not 
corrected for ellipticity. In addition, we tested a simple 
/?-model, with a shape matching that of the diffuse emis- 
sion, against the data. In both cases, we added an ad- 
ditional (constant) term to account for unrelated "back- 
ground" sources. We were able to obtain good fits to the 
data with both models (null hypothesis probabilities of 
56% and 28%, respectively), for which the inferred num- 
bers of background sources agreed within errors. Fig 
shows the data and the best-fit de Vaucouleurs model. 
For our de Vaucouleurs fit, we found a total of 130018qq° 
background sources per degree, equivalent to a total of 
1.9l]^'2 background sources in the D25 region. Taking 
our lowest measured flux to be the completeness limit, 
we would expect, (based up on the number of sources de- 
tected bv iTozzi et al.ll2001l in the Chandra Deep Field 
image), there to be between ~1400 and ~2600 sources 
per degree (for the "soft" and "hard" bands of these au- 



thors, corrected to our 0.3-7.0 keV band). Our data were 
clearly, therefore, in general agreement with the expected 
number of background sources from these observations. 

In addition, we accumulated a histogram of the az- 
imuthal source distribution. We binned the data into 128 
azimuthal bins, regrouping adjacent bins until there were 
at least 5 sources in each. We considered only sources 
within 2' of the galaxy centroid to minimize background- 
source contamination. Models were tested against the 
data using the Cash statistic, exactly as for the radially- 
binned data. The data clearly showed azimuthal varia- 
tion, as would be expected if the sources follow the op- 
tical light. Formally, we found that a model in which 
the source density is constant with angle gives a null 
hypothesis only of 1.8% and therefore was rejected. Al- 
ternatively, for a model in which the number of sources 
followed the optical light, the null hypothesis probability 
of 70% indicated excellent agreement with the data (see 
Fig©. 

4.8. Source completeness estimate 

A number of factors can influence the likelihood of de- 
tecting a given point source, such as source confusion in 
crowded regions, increased local background due to dif- 
fuse emission and the degradation of the Chandra point- 
spread off-axis. The importance of correcting the mea- 
sured XLF to account for incompleteness was demon- 
strated by Kim & Fabbiano (2003b), who showed that 
failure to account for it can lend erroneous significance 
to XLF breaks. In order to estimate the level of incom- 
pleteness in our source det ections we adopted a tech nique 
similar to that outlined bv lKim &: Fabbiaiiol l|2003bD . We 
ran a large number of Monte-Carlo simulations in which 
a small number of sources were added to the observed 
image of NGC 1332 in the 0.1-10 keV band and then 
measured the fraction of these sources which could be 
found with the detection algorithm. Since our adoption 
of additional energy-bands did not appreciably add to 
the number of detected sources, it was appropriate to 
consider only the 0.1-10 keV energy-band in this cal- 
culation. The added sources were all assumed to have 
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Fig. 6. — Left panel: The radial distribution of the point-sources, binned to ensure at least 5 sources per bin and shown with the best-fit 
model following the optical light. Right panel: The azimuthal distribution of sources, shown with the best-fit model, which also follows the 
optical light. 
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Fig. 7. — Representative plot of the detection probability for a 
point-source at a given luminosity as a function of measured counts. 
Also shown is our polynomial approximation to the function. 



spectra which were strongly peaked at 1 keV and appro- 
priate images were extracted from the Caldb PSF hyper- 
cube and added to the image (having been degraded with 
suitable Poisson noise). The sources were added with an 
assumed luminosity function, (dN/dS)o oc (although 
its adopted shape does not affect the results) and, in or- 
der to speed the computation, we added ~ 10 sources at a 
time to each image down to a limiting flux of 1 count per 
image. These were assumed to be distributed as the op- 
tical light (i.e. an elliptical de-Vaucouleurs profile with 
axis ratio and alignment matching that of the D25 el- 
lipse and with major-axis 0.55"; BC). Although ideally 
sources would be added one-at-a time so as not to intro- 
duce any biases, this expedient is not strictly necessary 
provided the incidence of overlap between added sources 
is low. After a total of 400 simulations, we were able 
to measure the function Pd(S), which is the probability 
that a source with exactly S counts is detected. Due to 
the finite number of test sources, we binned them into a 
number of '--^logarithmically-spaced luminosity-bins and 
computed the fraction of sources detected within each 
bin (Fig[7|). Pd(S) is given by the ratio of the measured 
and simulated differential XLF. We were then able to 



approximate Pd(S) by a fourth-order polynomial in the 
range S=l-50 counts (also shown in FigO. We consid- 
ered only simulated sources lying within the D25 ellipse, 
since we later computed the XLF only in this region. 

5. THE X-RAY LUMINOSITY FUNCTION 

In order to minimize the possible impact of background 
source contamination (a strong source of potential er- 
ror given field-to-field variations in the number of back- 
ground objects), when computing the XLF we consid- 
ered only those sources within the D25 region. We fit- 
ted the differential luminosity function, separately con- 
sidering the XLF derived from our two different flux- 
estimation techniques. Accordingly, the sources were 
binned as a function of luminosity in such a way that 
the bin size increased geometrically (progressively scal- 
ing by a factor 2.0). We experimented with different 
scaling factors but found our results were relatively in- 
sensitive to the adopted binning. Since our measurement 
included no upper luminosity cut-off, it is significant that 
we detected no sources with luminosities greater than 
~ 12 X 10'^^ erg s~^. Therefore, an additional data-bin 
containing zero sources and covering an extremely wide 
luminosity range (from Lx— 12-10000 xlO'^* erg s~^) 
was added to our data (the fit statistic is extremely in- 
sensitive to the precise upper bound of this bin, provided 
it is larger than ~ 100 x 10^^ erg s~^). Including this bin 
leads to steeper XLF fits than omitting it. To fit the XLF 
we minimized the Cash statistic since our data were dom- 
inated by Poisson errors. The fit quality of the model was 
subsequently determine via Monte-Carlo simulations. To 
take account of flux measurement errors, we repeated the 
fit 100 times for each model, having first added appro- 
priate Gaussian noise to the measured Lx of each source. 
If any resulting Lx would be below our lowest measured 
luminosity (assumed to be the source detection limit), 
it was recomputed. In order to test the reliability of 
this fitting procedure, we simulated 100 luminosity func- 
tions from our best-fit model and attempted to recover 
the best-fit parameters. We found good agreement be- 
tween the inferred and expected results and confirmed 
the expected statistical scatter, giving us confidence in 
our ability to determine parameters and estimate error- 
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Fig. 8. — The differential luminosity-function of the data. Error-bars, estimated from the Gehrels' approximation, have been included 
on each point only as a guide to the eye since the Cash statistic was used in fiting. Also shown are the best-fit completeness-corrected 
power law model (solid line) and with the best-fit uncorrected simple power law model (dashed line). This demonstrates the impact of 
incompleteness on measuring the XLF. 



bars reliably, provided the simulations were representa- 
tive, which we believe. 

We tested a number of models against the data, in 
each case adding a separate component to account for 
the point-source population due to background sources. 
This term had a power law shape, i.e. 



dL V1038erg s-iy ' ' 

We fixed the slope, /?, and the amplitude, K, to val- 
ues co nsistent with the hard-band sources of Tozzi et al.l 
l)200Hl . i.e. 1.92 and 2.1 sources per 10^^ erg s~^, re- 
spectively (correcting fluxes from the 2.0-10.0 keV band 
used by these authors to our 0.3-7.0 keV band). This 
term did not contribute very significantly to our fit and 
corresponds to ~ 3.5 sources in total in the D25 re- 
gion. We found that allowing the normalization to be 
free, or adopting the sha pe and no rmalization appro- 
priate for the soft band of iTozzi et a l. (2001) made lit- 
tle difference to our overall fit. Without applying any 
form of completeness correction to our data, we found 
that a single power law fit (with fluxes estimated from 
spectral-fitting) gave a very poor fit to the data (having 
a null hypothesis probability of 0.3%). Adding a break in 
the model at Lbrcak = 2.3i|!)'7 x 10^^ erg s^^substantially 
improved the fit quality (giving a null hypothesis prob- 
ability of 42%), although the power law slopes were 



poorly-constrained (e.g. (3 = 2.61q y for Lx> Lbreak)- To 
correct for source detection incompleteness we modified 
the simple power law by multiplication with our com- 
pleteness estimate function (i.e. we corrected the model 
rather than the data, to allow us to continue using Pois- 
son statistics). We obtained an excellent fit with the cor- 
rected single power law model (the null hypothesis prob- 
ability now being 40%). The best-fit parameters were 
(3 = 2.3 ± 0.4 and amplitude, K=60 ± 25, and no break 
was required. This indicates that incompleteness was re- 
sponsible for the break measured in the uncorrected data. 
No further improvement in the fit could be obtained by 
adding a break to the power law model, and the resulting 
fit-parameters were very poorly-constrained. The mea- 
sured value of was in good agree ment with that ob- 
tained by Kim fc Fabbianol l)2003bD for a sample of 14 
early- type galaxies. However, these authors computed 
source fluxes based on a counts-to-flux conversion, so this 
is not strictly a fair comparison. 

We separately considered the XLF where fluxes were 
obtained from the count-rate, rather than the spectra. 
This allowed us to compare our results more directly with 
other authors, who tend to adopt this procedure. Ex- 
actly as above, we found that a broken power law gave a 
much better fit than a simple power law to the data, but 
this effect vanished when proper account is made for the 
source detection incompleteness. In this case, however. 
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we obtained a somewhat steeper slope for our luminosity- 
function {f3 — 2.7 ±0.5; K was correspondingly altered to 
87 ± 33 ), although the two f3 values agree within errors. 
The differential XLF and the best-fit model are shown in 
Fig. IHl along with the uncorrected single power law fit. 
It is not surprising that the XLF derived from fitting the 
spectra was less steep, since the principal effect of the fit- 
ting was to uncover sources with more than Galactic Nh, 
thereby populating more the higher luminosities. Since 
some of these "absorbed" sources may have been back- 
ground AGN and since the background XLF we adopted 
were computed from counts-to-flux conversioii, it is prob- 
ably more correct to adopt our f3 = 2.7 ± 0.5 result. In- 
triguingly, this result is marginall y (~ 2.5a) steeper than 
the 13 = 2.0 value found by Kim & Fabbiano ( 2G03b) but 
is in good agreement with their higher-luminosity slope 
when they fitted a broken power law to the combined, 
incompleteness-corrected XLF. To verify that our slope 
was not biased by our use of data during the mild flaring, 
we recomputed all fluxes having excised data during the 
flares. The slope of the XLF was unchanged. 

As a final note, we also attempted to fit the XLF using 
two alternative fitting methods: a max imum-likelihood 
meth od appropriate for unbinned data ijCrawford et all 
|1970|) and minimization of the Kolmogorov-Smirnov (K- 
S) test statistic. In both cases we obtained less strin- 
gent constraints than for our binned-data analysis above. 
Since, by its nature, the K-S test is relatively insensitive 
at the sparsely-populated highest and lowest luminosi- 
ties, that statistic in particular was much less reliable at 
determining the quality of the fit. 

6. DISCUSSION 

We have found a total of 73 point-sources within the 
Chandra S3 field of NGC 1332, including 37 which lie 
within the D25 isophote and none of which is coincident 
with the galaxy centroid. We expect a total of 1.9l]^'2 
unrelated background or foreground sources to be within 
D25, so that the vast majority of these sources are as- 
sociated with the galaxy. The spatial distribution of 
the sources was consistent with the optical isophotes, as 
may be expected from a population of sources of stellar 
origin w hich has not been recently disturbed by merger 
activity l|Zezas et al.l 120031) . In total 30% of the X-ray 
sources were associated with globular clusters, and 4.4% 
of GCs contained an X-ray source, both of which results 
are consistent with other early-type galaxies. The phe- 
nomenology of the individual sources was consistent with 
expectations for an LMXB population. Two sources had 
soft spectra, consistent with high-state black-hole bina- 
ries and two sources were variable during the observation. 
Intriguingly one source was extended, although its radial 
profile was consistent with two unresolved point-sources. 

Without taking into account the impact of source de- 
tection incompleteness at the lowest luminosities (where 
source confusion and increased background due to dif- 
fuse gas reduce the probability of detecting a source), 
we found that the XLF of the point-sources within 
D25 can be fitted by a power law with a break at 
Lx~ 2 X 10'^^ erg s~^. Similar results have been re- 
ported in other earl y-type galaxies (e.g. iSarazin et al.l 
l2001t iBlanton et a l. 2001). However, when we made 
corrections for the expected level of incompleteness, we 
found that the data were extremely well-fitted by a sin- 



gle power law, with a differential slope, f3 = 2.7 ± 0.5. 
Assuming that the XLF remains unb roken down to 10^ ^ 
erg s~^ (as seen for sources in M31: IKong et alJl200l . 
we estimate a total luminosity of 180 ± 70 x 10'^^ erg s~^ 
due to unresolved point-sources in the D25 region. Com- 
paring this to the B-band luminosity of the galaxy, ob- 
tained from the face-on, extinction-corrected B-band lu- 
minosity listed in the LEDA'^ catalogue, we obtained 
Llmxb/Lb = 0.8 ± 0.3 X 10^° erg s'^LgQ. In addition, 
we compared with the K-band luminosity (adopting 
the K20 magnitude obtained by 2MASS), and obtained 
Llmxb/Lk = 0.13 ± 0.05 x lO^" erg s-'^L-^^. Both of 
these values were in excellent agreement with obser- 
vations of other early-type galaxies (Kim & Fabbianol 
i2003b). 

The slope of our XLF fit closely resembles that found at 
high luminosit i es in a variety of other early-type galaxies 
l|g;arazin et a,].ll200Tt IB]a,nto;77nni200Tt rColbert et alJ 

l2003f) . A slope of /3 ~ 2.7 is also strikingly similar to ob- 
servations of a primarily old stellar population in M31, 
but which does n ot extend to such high luminosities 
l)Kong et al.ll200.^ . The slope is, however, marginally 
(~ 2.5cr) steeper than the incompleteness-corrected sin- 
gle power law fits found in a recen t survey of early-type 
galaxies ijKim fc Fabbianol l2003b(l . Intriguingly, these 
authors reported some evidence of a break at 5 x 10^^ 
erg s~^, above which /3 ~ 2.7 is ob served. A similar re- 
sult has also been found for M84 ("Finog uenov fc JonesI 
f2002). A break at this luminosity should be observable 
within our data, and its apparent absence is intriguing. 
Whilst it is certainly true that statistically acceptable 
single po wer law fits were found for most galaxies in the 
sample of iKim fc Fabbialiol l)2003bf) . the resulting slope 
{(3 ~ 2.0) may simply be an "average" of the slopes above 
and below the break. The steep slope seen in NGC 1332 
might therefore suggest that such a break is absent in this 
galaxy. It must be stressed, however, that there are some 
differences in the analysis procedures adopted by these 
authors and those used in the present work, which may 
simply weight our fit somewhat more heavily towards 
the high luminosity sources (especially our inclusion of 
a high- luminosity, unoccupied bin). Due to the small 
numbers of sources, however, the error-bars on the best- 
fit slope in NGC 1332 are large, so that the disagreement 
is only at ~ 2.5f7. It remains to be seen whether com- 
parable results are observed in any other galaxies con- 
taining significantly more X-ray point sources. Nonethe- 
less, if the absence of the break in our XLF is confirmed, 
it is an interesting result, which casts further doubt on 
there being a "universal" break in early-type galaxies 
corresponding to the division between neutron-star and 
black- hole binary systems. 

The origin of the XLF break is still somewhat un- 
clear. Although it was originally suggested that a break 
might mark the division between neutron-star and black- 
hole binary systems (jSarazin et al... 2001). the dispar- 
ity between the reported break luminosities in several 
galaxies and the Eddington limit for a 1.4M0 (canoni- 
cal) neutron-star suggests another origin, perhaps in the 
star- formation history of the galaxy (e.g. iJeltema et all 
12003: Wu 2001). In that case, however, we might expect 
to see a strong correlation between the age of the sys- 

'* |http://Ieda.univ-Iyonl.fr/| 
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tern (since the last major me rger) and the XLF shape ; 
comparing the slopes found bv lKim fc Fabbianol l)2003hD . 
including that for NGC 7 20, and galaxy ages given by 
iTerlevich fc Forbe^ l)2002f) we find no evidence of a cor- 
relation. However, this includes only 8 sources and will 
need to be confirmed with a larger sample. 

We detected two ULX within the D25 region of 
NGC 1332, although neither of these sources have Lx> 
2 X 10'^^ erg s~^. Such numbers are consistent with other 
galaxies of similar morphology (Humohrev et al. 200!j). 
The unbroken XLF of late- type galaxies suggests that the 
ULX associated with star-formation can be attributed 
to the high-luminosity tail of the HM XB XLF, w hich 
is flatter than the LMXB XLF (Grimm eral]l2003|) . It 
may be that, in a similar way, the ULX in early- type 
galax ies simply represent the high-Lx tail of the LMXB 
XLF. IIrwin et al.l l)2003bD have argued that sources with 
Lx> 2 X 10'^^ erg s~^ are extremely rare in early-type 
galaxies, suggesting that the most luminous ULX may 
therefore represent a population almost entirely alien 
to such galaxies. However luminous ULX have been 
reported in NGC 1399 and NGC 720 which are associ- 
ated with GCs (hen ce u nlikely to be backgr ound ob- 
jects; lA^ielinreraOllMIl EsEema^^ so they 
are clearly not entirely absent and therefore it is of in- 
terest to determine if their paucity simply reflects the 
steepness of the XLF. In NGC 1332, since our fitting took 
account of the absence of sources at high luminosities, 
the absence of very luminous ULX is consistent with the 
high-luminosity tail of the XLF. However, we also com- 
parcd the ratio of the number of sources measured by 
llrwin et all l)2003bD in the 1-2 xlO^^ erg s'^ and the 



> 2 X 10'^^ band, to that expected assuming our best-fit 
XLF slope iP = 2.7). We found a likehhood of - 1% 
of reproducing the data from the model, although al- 
lowing P to steepen within its error-bars increased the 
likelihood to ~ 5%. Nonetheless this suggests the possi- 
bility of a change in the XLF at 1-2x10 ^^ erg s~^. In- 
triguin gly, an XLF break was reported bv lJeltema et alJ 
l)2003j) around this luminosity. 
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TABLE 1 

Sources detected 



Src 


Name 


Rate 
(10-4 s-i) 


(lO^^erg s-i) 


^rate 
(lO^^erg s-i) 


AR 
(') 


In D25 


Comment 


1 


CXOU J032636. 5-211809 


81. 


58. ± 5. 


60. it 5. 


4.9 


no 


R 


2 


CXOU J032614.8-212126 


95. 


53. ± 4. 


46. it 3. 


1.4 


no 


A, R, NG 


3 


CXOU J032604.4-212258 


45. 


36. ± 5. 


24. it 3. 


4.1 


no 


A 


4 


CXOU J032620.7-212151 


11. 


23.+I- 


5.3 it 1.3 


1.9 


no 


A, NG 


5 


CXOU J032632.8-212306 


6.0 


— 3. 


8.6 it 3.6 


4.7 


no 


A, C 


6 


CXOU J032618.1-212014 


18. 


12.+3- 

— 2. 


12. it 2. 


0.23 


yes 


A, B, R, NG 


7 


CXOU J032601. 7-211726 


21. 


12. ± 2. 


12. it 2. 


4.5 


no 




8 


CXOU J032616.7-211954 


19. 


12. ± 2. 


11. ± 2. 


0.26 


yes 


NG 


9 


CXOU J032618. 7-212031 


17. 


7.8 ± 1.8 


8.4 ± 1.6 


0.51 


yes 


NG 


10 


CXOU J032606. 8-212016 


12. 


■7 0-1-2.0 


6.1 ± 1.4 


2.4 


no 




1 1 
i 1 


*^jV^ U J UoZDiO.y-Z iZUUD 


0.0 


0. / it Z.D 


0.0 ± Z.o 


U.uy4 


yes 






t^jjVV^ U J UoZOiZ.O-Z iZo4Z 


D.Z 




Q -1- 1 
0.0 It l.Z 


0. / 


no 


A 




v_>jVW u J uoZDUy • ( -z izzi i 


1 


4- 1 

o.y It i.y 


4.D It 1.0 


V 
Z. / 


no 




1 A 


r^VOTT TflQO^^I /I fi 01100*3 
(^AW U J UoZD14.D-Zi lyUo 


1 n 
iU. 


0.0 It l.D 


^ 1 "3 
D.Z ± 1.0 


1 
l.Z 


no 




15 


CXOU J032d15.z-z1z001 


9.3 




5.6 it 1.5 


0.50 


yes 




16 


CXOU J032622. 6-212418 


8.5 


_ p-t-2.0 
0.t)_j 2 


4.7 it 1.4 


4.3 


no 




17 


CXOU J032615. 2-211924 


2.9 


A n-l-2.1 


1.5 it 0.8 


0.87 


yes 


A 


18 


CXOU J032617.1-212011 


2.8 


,1 tr-l-2.1 

4-5Ii.6 


1.9 it 2.2 


0.082 


yes 


G 


19 


CXOU J032617.6-212013 


2.8 


^•^-1.6 


2.0 it 2.4 


0.13 


yes 


G 


20 


CXOU J032622.8-211917 


7.3 


4 4+1-8 


3.7 it 1.1 


1.5 


no 




21 


CXOU J032603. 1-21 1921 


7.4 


4 3+2-3 


3.7 it 1.1 


3.4 


no 


A, B 


22 


CXOU J032631. 3-212339 


5.8 


3.9 it 1.8 


3.2 it 1.2 


4.8 


no 




23 


CXOU J032620. 0-212023 


5.7 


3.4 it 1.3 


3.0 it 1.1 


0.69 


yes 


NG 


Z4 


t^AWU JUOZDID.O-Zl ±yoo 


A 7 


It 1.0 


z. / It 1.4 


U.Z / 


yes 




25 


CXOU 1032622 4-212132 


7.1 


Q 9-H-6 

•^•^-0.8 


3.5 it 1.1 


1.8 




NG 


26 


CXOU T03261 3 1-21 2020 


4.7 




2.3 it 0.9 


0.99 


yes 


NG 


27 


CXOU 1032618 7-211814 


5.9 


3 0+1'^ 


^ I -J- 11 


1.9 




Y 


28 


CXOU 1032605 7-212058 


5.3 


2 7+1-1 
^•'-0.9 


K + 1 no 


2.8 






29 


CXOU J032619.0-212030 


4.3 


2.4 it 0.9 


2.2 ± 0.9 


0.54 


yes 


G 


30 


CXOU J032619.4-212021 


3.8 


2.4 it 1.0 


2.0 ± 1.0 


0.53 


yes 


G 


31 


CXOU J032619.3-212052 


3.6 


2.3+?-I 


1.8 it 0.8 


0.88 


yes 


NG 


32 


CXOU J032619.9-212019 


3.0 


2 3+11 


1.9 it 1.00 


0.64 


yes 


NG 


33 


CXOU J032616.4-211953 


4.4 


9 0-H.2 


2.2 it 1.0 


0.31 


yes 


G 


34 


CXOU J032601.6-212021 


3.6 


2 3+11 


1.8 it 0.9 


3.6 


no 




35 


CXOU ,1032625.5-211622 


2.1 


9 Q-l-1-5 


1.8 it 1.1 


4.2 


no 


C 


36 


CXOU J032618.3-212011 


3.9 


9 o-H-5 


2.4 it 1.2 


0.25 


yes 


NG 


37 


CXOU J032621.6-211930 


3.5 


9 9-H0.8 


1.7 it 0.8 


1.2 


no 


NG 


38 


CXOU J032627.4-212235 


5.7 


2 2+^-4 


2.9 it 1.2 


3.4 


no 




39 


CXOU J032616.2-211953 


4.3 


2 1+^-2 
^•1-0.7 


2.2 it 1.0 


0.34 


yes 




40 


CXOU J032615.5-211941 


4.8 


2 1+1-1 


2.4 ± 1.0 


0.60 


yes 




41 


CXOU J032602.2-212140 


3.7 


2 0+1-3 


2.0 it 1.2 


3.8 


no 




42 


CXOU J032620.6-211713 


3.6 


1 9+1-3 
l-^-0.5 


1.8 it 0.8 


3.0 


no 




43 


CXOU J032613.9-211944 


3.6 


1 S+1-4 


1.8 it 0.8 


0.87 


yes 




44 


CXOU J032610.3-211842 


3.1 




1.6 it 0.8 


2.1 


no 




45 


CXOU J032620.1-211958 


2.3 


1 7-1-0.9 
l-'-0.7 


1.2 it 0.7 


0.67 


yes 


NG 


46 


CXOU J032620.4-212112 


3.0 


1 6+°-^ 
l-"-0.7 


1.5 it 0.8 


1.3 


yes 


NG 


47 


CXOU J032631.3-212210 


3.2 


1 6+1-1 


1.6 it 0.9 


3.8 


no 




48 


CXOU J032624.2-212040 


2.0 


1 6+°-'^ 


1.2 it 0.8 


1.7 


yes 




49 


CXOU J032620.6-212317 


2.2 


1 5+1-° 


1.2 it 0.8 


3.2 


no 




50 


CXOU ,1032630.9-211824 


3.3 


1 ,;-|-0.9 
1-^-0.6 


1.7it0.8 


3.6 


no 




51 


CXOU J032606.7-212139 


3.0 


1 4+1-2 


1.5 it 0.8 


2.9 


no 


V 


52 


CXOU J032621.9-212023 


3.3 




2.1 it 1.00 


1.1 


yes 


NG 


53 


CXOU J032614.6-211955 


1.7 


1 o-Hl.O 

i-^-o.s 


0.93 it 0.70 


0.65 


yes 
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TABLE 1 
Sources detected (contd.) 



Src 




Name 


Rate 
(10-*s-i) 


J spec 

(10^** erg s~i) 


T rate 

(lO^^erg s-l) 


AR 

w 


In D25 


Comment 


54 


CXOU 


J032620. 5-212041 


2.5 


1 O+0.7 
^•"^-0.6 


1.3 lb 0.7 


0.94 


VPS 


G 


55 


cxou 




2 1 


-'-•^-0.4 


1 1 i 7 


0.79 


yes 




56 


CXOU 


J 032625. 1-21 2043 


2.0 


1 1+1-0 
1-1-0.4 


1.1 it 0.7 


1.9 


yes 




57 


cxou 


J 032621. 8-21 1930 


2.6 


1 1+0.5 
1-1-0.7 


1.3 it 0.7 


1.2 


no 


NO 


58 


cxou 


J032609. 2-212156 


1.5 


1 1+0.5 
1-1-0.7 


0.81 it 0.66 


2.6 


no 




59 


cxou 


J032622. 8-21 1620 


1.4 


1 1+0.8 
1-1-0.5 


0.76 it 0.64 


4.0 


no 




60 


cxou 


J032623. 1-21 1945 


1.6 


1 1+0.7 
1-1-0.5 


0.79 it 0.60 


1.4 


no 


NG 






TIT^Qfilfi ^ 911090 


1 7 


1 n -t- n 


u.oo in u.uo 


U.O-L 


yes 




62 


cxou 


J032634. 7-211755 


2.0 


96+l L 


1.1 it 0.7 


4.6 


no 




63 


cxou 


J032615. 4-211948 


2.1 


96+l L 
"-^"—0.34 


1.0 it 0.7 


0.53 


VPS 




64 


cxou 


J032612.8-211960 


2.1 


n 93+1'^ 
u.»o_g 33 


1.1 it 0.7 


1.0 


yes 




65 


cxou 


J032615. 7-212038 


2.6 


0.93 it 0.57 


1.3 it 0.8 


0.63 


yes 


NO 


66 


cxou 


J032620.8-211957 


1.4 


0-88tr4^ 


0.73 it 0.61 


0.83 


yes 


G 


67 


cxou 


J032621. 5-212029 


1.3 


0-821^9 


0.68 it 0.62 


1.1 


yes 


G 


68 


cxou 


J032616.6-211927 


1.6 




0.79 it 0.65 


0.69 


yes 




69 


cxou 


J032604.7-211613 


1.7 


0.781^4^ 


0.89 it 0.66 


4.8 


no 




70 


cxou 


J032607.9-211828 


1.5 


0.751^4^ 


0.77 it 0.61 


2.7 


no 




71 


cxou 


J032618.8-212148 


1.4 


Q co+0.56 

U-DO_g 3g 


0.69 it 0.60 


1.7 


no 


R(?), NG 


72 


cxou 


J032614.2-211950 


1.4 


Qgo+0.57 


0.75 it 0.67 


0.76 


yes 




73 


cxou 


J032615.7-212017 


1.2 




0.61 ±0.59 


0.39 


yes 


NG 



Note. — Source list for all point-sources detected on the ACIS S3 chip, excluding the source coincident 
with NGC 1331. Two luminosity estimates are given, L^'^'^ and L^''^ which are, respectively, the luminosity 
estimated from the spectra and from the count-rate (exposure-corrected). AR is the distance of the source 
from the galaxy centroid (in arcmin). The count-rate column quotes count-rates not corrected for exposure- 
map variations. The comment abbreviations are as follows: A — absorbed spectrum, B — soft spectrum, C — 
close to chip edge (flux may be a unreliable), E — extended, G — coincident with a globular cluster, NG — in 
the WFPC2 fov but not coincident with a GC, R — also seen with Rosat and V — variable. 



TABLE 2 
Composite source spectra 



Model X^/dof Nh T kT 

10^2 cm- 2 (keV) 



All sources 


power law 122/119 
115.3/118 


0.0223 
0.07 it 0.03 


1.40 it 0.05 
1.5 ± 0.1 


15t^ 


bremsstrahlung 115/119 


0.0223 




114/118 


0.03 it 0.02 




12tl 


Sources in D25 only 


power law 50.0/46 
50.0/45 


0.0223 
< 0.075 


1.55 it 0.10 
1.6 ± 0.1 




bremsstrahlung 50.8/46 


0.0223 




Silt 


49.2/45 


< 0.031 




10±3 



Note. — Results of fitting the composite spectra of 
all the point sources on the ACIS-S3 chip and, separately, 
only those sources in the D25 region. Njj values shown in- 
clude the Galactic column density, which was 0.0223xl022cm~2 
IDickev fc Lockman|[T99Cl) . Where no error-bars are shown the 
parameter was fixed. Errors are 90% confidence intervals. 



